arXiv: 1503.05565vl [cond-mat.quant-gas] 18 Mar 2015 


Measuring the dynamic structure factor of a quantum gas undergoing a structural 

phase transition 

Renate Landig,^ Ferdinand Brennecke,^ Rafael Mottl/ Tobias Donner*,^ and Tilman Esslinger^ 

^Institute for Quantum Electronics, ETH Zurich, CH-8093 Zurich, Switzerland 
^ Physikalisches Institut, University of Bonn, Wegelerstrasse 8, 53115 Bonn, Germany 

(Dated: March 20, 2015) 

The dynamic structure factor is a central quantity describing the physics of quantum many- 
body systems, capturing structure and collective excitations of a material. In condensed matter, 
it can be measured via inelastic neutron scattering, which is an energy-resolving probe for the 
density fluctuations. In ultracold atoms, a similar approach could so far not be applied due to the 
diluteness of the system. Here, we report on a direct, real-time and non-destructive measurement of 
the dynamic structure factor of a quantum gas exhibiting cavity-mediated long-range interactions. 

The technique relies on inelastic scattering of photons, stimulated by the enhanced vacuum field 
inside a high finesse optical cavity. We extract the density fluctuations, their energy and lifetime 
while the system undergoes a structural phase transition. We observe an occupation of the relevant 
quasi-particle mode on the level of a few excitations, and provide a theoretical description of this 
dissipative quantum many-body system. 


An interacting quantum many-body system can be 
characterized by analyzing its response to a weak per¬ 
turbation. In the framework of linear response theory a 
key quantity is the dynamic structure factor, which is the 
Fourier transform of the spatial and temporal density- 
density correlations mm- Knowledge of the dynamic 
structure factor provides a complete picture of the emerg¬ 
ing quasi-particle modes [3], their excitation energy, life¬ 
time and mean occupation number. These quasi-particle 
modes determine the collective density fluctuations of the 
system and may also characterize the critical behavior in 
the vicinity of a phase transition |3]. For example, in a 
long-range interacting system a structural phase transi¬ 
tion can be driven by a roton-like mode softening I3ISH3, 
which is expected to show up as a thermally enhanced 
peak in the dynamic structure factor [9]. 

In solid state systems, the dynamic structure factor 
S'(k, w) can be measured by illuminating a sample with 
a beam of neutrons m or x-rays m, and analyzing 
the inelastically scattered particles with regard to their 
change in energy Huj and momentum fik. In dilute 
quantum gases, a measurable signal has only been ob¬ 
tained from photons elastically scattered off a density- 
modulated sample [MIS]. Direct detection of inelasti¬ 
cally scattered photons into free space [T6| , in analogy to 
neutron scattering, is however hindered by a vanishingly 
small signal m, see Fig. lA. 

A technique measuring the spectral response function 
of a quantum gas, i.e. the dynamic structure factor at 
zero temperature [18] , is Bragg spectroscopy [H] [20] . It 
is based on stimulated rather than spontaneous inelastic 
scattering of photons between two laser beams. There¬ 
fore the transfer of momentum and energy to the atomic 
cloud is predetermined by the angle and frequency differ¬ 
ence between the beams, and is typically measured via 
destructive absorption imaging. A complementary detec¬ 
tion method analyzes the change in light field intensity in 
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Figure 1. Scheme for measuring the dynamic struc¬ 
ture factor in a qnantum gas. (a) An incident laser beam 
(red) with wave vector ki is spontaneously scattered at an 
atomic cloud (blue) into a free-space mode with wave vector 
kf. Analyzing the scattered photons as a function of their fre¬ 
quency shift Lj{ — LJi and magnitude yields the dynamic struc¬ 
ture factor. The resulting signal for dilute quantum gases is 
vanishingly small, as indicated by the small size of the ar¬ 
row pointing towards the detector, (b) Atoms placed into an 
optical high-finesse resonator (light blue) feel a strongly en¬ 
hanced vacuum field (golden). Their spontaneous scattering 
rate into the mode kf can hereby be increased by orders of 
magnitude, resulting in a detectable signal for the dynamic 
structure factor. 


one of the two Bragg beams |21j . and could in principle 
be extended with the help of cavities to be only weakly 
perturbative [52]. However, all these methods measure 
the linear response of the gas upon a perturbation and 
are insensitive to thermally excited quasi-particles. A 
different approach, in-situ imaging, has been used to ex¬ 
tract the temperature-dependent static structure factor 
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S'(k) = f dLjS(k,cj) of a two-dimensional gas [231. 
ilar to the analysis of noise correlations from images of 
ballistically expanded ultracold gases |2l], this approach 
gives no access to the quasi-particle spectrum, i.e. the 
temporal dynamics. 

Here, we present a non-destructive, direct measure¬ 
ment of the dynamic structure factor at distinct k-vectors 
in a Bose-Einstein condensate (BEC) undergoing a struc¬ 
tural phase transition induced by cavity-mediated long- 
range interactions. We place a BEC into an ultra-high 
finesse optical cavity [231123] and illuminate the atoms 
with a transverse laser field. The enhanced vacuum field 
inside the optical resonator [23 increases the spontaneous 
inelastic scattering of photons into the cavity mode by 
several orders of magnitude, so that photons leaking out 
of the cavity mode give rise to a detectable signal and 
access to the density correlations in real time, see Fig. 
IB. Therefore density fluctuations in the gas are mapped 
onto fluctuations of the light field, which then can be 
directly accessed. 


RESULTS 



System description 

In our experiment, the transverse laser field acts si¬ 
multaneously as a pump field controlling the long-range 
interactions This field has a frequency Wp and a 

wavevector kp, and is in a standing-wave configuration 
directed perpendicularly to the cavity mode. It is far de¬ 
tuned from atomic resonance to avoid electronic excita¬ 
tion of the atoms. At the same time, it is detuned by only 
a few cavity linewidths from the cavity resonance, which 
enables vacuum-stimulated scattering of pump photons 
into the cavity mode at wavevector kc. These two-photon 
processes mediate long-range atom-atom interactions in 
the BEC, giving rise to a roton-like mode softening |3] 
and a structural phase transition mm- The same two- 
photon processes are exploited for detection. The light 
scattered from the transverse pump field into the cavity 
mode can be regarded as a superposition of all field am¬ 
plitudes scattered by the individual atoms. It thus carries 
information on the density-density correlations of the gas 
at the wavevectors kcb = ±kp±kc, which are determined 
by the underlying two-photon processes [2]. Within the 
cavity linewidth, which is two orders of magnitude larger 
than the frequency of the relevant quasi-particle exci¬ 
tation, the energies of the photons stimulated into the 
vacuum mode are not fixed. This is in contrast to Bragg 
spectroscopy, where the energy of the quasi-particle exci¬ 
tations that are created during probing is determined by 
the frequency difference of the classical fields driving the 
two-photon processes. The spectral analysis of the light 
field leaking out of the cavity thus gives us direct access 
to the dynamic structure factor at finite temperatures. 


Figure 2. Power spectral density and dynamic struc¬ 
ture factor, (a) The power spectral density PSD of the 
light field leaking out of the cavity is shown as a function of 
frequency shift oj with respect to cdp and relative transverse 
pump power P/Pct (average over 147 experimental repeti¬ 
tions). Two sidebands are visible, corresponding to the in¬ 
coherent creation (ca < 0) and annihilation (tj > 0) of quasi¬ 
particles. The energy of these quasi-particles vanishes towards 
the critical point. At the phase transition, a strong coherent 
field at the pump frequency appears (a; = 0). We attribute 
the broadened feature around a; = 0 to residual low-frequency 
technical noise in our system. Note the large dynamic range 
of the data on the logarithmic scale. The panels (b-d) show 
the normalized dynamic structure factor ^(kcb,^;) for three 
different values of P/Pcr (see dashed lines in upper panel), 
derived from the power spectral density PSD{uj). While the 
position and width of the sidebands give direct access to the 
energy and lifetime of the quasi-particles, the sideband asym¬ 
metry can be used to determine the occupation of the quasi¬ 
particle mode. Red line shows a fit to the sidebands with our 
theoretical model to extract these properties. 


As previously described [21 121 IHl [20] , we trap a BEC 
of = 1.0(1) X 10® ®^Rb atoms at the center of 
an ultrahigh-finesse optical Fabry-Perot cavity and il¬ 
luminate it by the transverse pump field. The cavity- 
mediated interaction leads to the formation of a quasi¬ 
particle mode, which is a superposition of the collective 
momentum excitation at wavevectors kcb of the BEC and 
a tiny admixture of photons inside the cavity (see Sup¬ 
plementary Note 1). Neglecting atom-atom collisions, its 
energy Hujs equals in the limit of zero pump power P 
the bare energy hujQ = {2m) of a single momen¬ 

tum excitation and decreases with increasing power |2], 
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Figure 3. Density fluctuations and density modulation. 

The variance of density fluctuations (filled blue symbols) and 
the square of the density modulation (open black symbols) of 
the long-range interacting quantum gas is shown as a func¬ 
tion of relative pump power P/P^r- The fluctuation data is 
extracted from the dynamic structure factor (see Fig. 2) by 
integrating over the fit function describing the sidebands and 
is proportional to the static structure factor. The coherent 
density modulation is calculated from the power spectral den¬ 
sity at the zero frequency bin. The vertical error bars display 
the statistical error (s. d.) from the fit, while the horizontal 
error bars display the standard deviation in our determination 
of the critical point. The inset displays a double logarithmic 
plot to demonstrate the scaling behavior of the variance of the 
density fluctuations against the distance to the critical point, 
expressed as the Hamiltonian coupling parameter A (see Sup¬ 
plementary Note 1). From a linear fit, we find critical expo¬ 
nents of 0.7(1) and 1.1(1) on the normal and self-organized 
side, respectively. The open symbols in the inset are used for 
the fitting. 

where m denotes the atomic mass. Due to this mode 
softening, the energy Hujs of the quasi-particle mode ap¬ 
proaches zero at a critical pump power Per, which leads to 
a phase transition from a normal state with a flat den¬ 
sity distribution to a self-organized state with checker¬ 
board density modulation. The emergent density struc¬ 
ture leads to elastic scattering of transverse pump light 
and a macroscopic population of the cavity mode. An 
order parameter of this phase transition is the expecta¬ 
tion value of the operator 0 describing the overlap of the 
atomic density p{r, t) and the checkerboard mode struc¬ 
ture, 0(t) = / d^r/5(r, t) cos(kcr) cos(kpr). 

Dynamic structure factor 

As the cavity decay rate k = 27r x 1.25 MHz is more 
than two orders of magnitude faster than the evolution 
rate Wg of the coupled system, the light field d{t) in¬ 
side the cavity adiabatically follows the order parameter, 
a{t) oc 0(t) [2]. The frequency spectrum of the light leak¬ 
ing out of the cavity thus reveals the temporal and spa¬ 
tial Fourier transform of the atomic density correlations, 
evaluated at one particular wavevector. Specifically, the 


dynamic structure factor of the system at wavevector keb 
is related to the cavity field according to 

5(k,b,a;) = ’ 

( 1 ) 

where w is the frequency shift of the cavity output field 
from the pump light frequency Wp due to inelastic scat¬ 
tering. PSD{uj) is the power spectral density of the 
intra-cavity light field with mean coherent field ampli¬ 
tude a = (d), and rj is the two-photon Rabi-frequency 
of the scattering process, proportional to Vp. Ac/27r = 
15.1(2) MHz is the detuning between the pump laser fre¬ 
quency and the dispersively shifted cavity resonance (see 
Supplementary Note 1). 

To analyze the light field leaking out of the cavity, we 
use a balanced heterodyne detection scheme [SU]. Figure 
2 shows the power spectral density PSD{uj) of the light 
field as we linearly increase the transverse pump power 
P across the critical point. The rate of change of P/Per 
is a few Hertz, such that the system can be assumed 
to adiabatically follow its steady state throughout the 
measurement [2 [30]. The small panels in Fig. 2 show 
examples of ^(keb,^) for different values of P/Pa , con¬ 
verted via Eq. ([H and normalized to unity for the non¬ 
interacting case [3T]. The data reflects the microscopic 
processes taking place: Pump photons of frequency Wp 
inelastically scattered at the atomic ensemble will be 
shifted in their frequency. They become visible as red 
(blue) sideband at frequency Wp — Wg (wp + Wg) if they 
create (annihilate) a quasi-particle. We observe the cor¬ 
responding sidebands whose frequency shift tends to zero 
when approaching the critical point at P/Per = 1 from 
either side of the phase transition. These density fluc¬ 
tuations can be distinguished from a checkerboard den¬ 
sity modulation at keb at which pump photons will be 
elastically scattered without a frequency shift, visible at 
w = 0. Intuitively this light field arises from scattering 
at Bragg planes in the density modulated cloud. The 
transverse pump power P influences not only the effec¬ 
tive long-range interactions in the system, but also the 
measurement process itself. For increasing pump power, 
the measurement imprecision due to the shot noise of 
the transverse pump field becomes less relevant, as can 
be seen from the decreasing background level of S'(kcb, w) 
|32) . The influence of the inevitable measurement back- 
action will be discussed further below. 

The total power of elastically scattered light is pro¬ 
portional to the square of the density modulation (0)^, 
and is displayed as open symbols in Fig. 3. This co¬ 
herent density modulation increases over five orders of 
magnitude while crossing the critical point. In the nor¬ 
mal phase {P/Per < 1), we also observe a weak field at 
the pump laser frequency (i.e. at a; = 0), which orig¬ 
inates from a small symmetry breaking held caused by 
the finite size of the system and residual scattering of 
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the transverse pump beam at the cavity mirrors [5J I30j . 

The total power of frequency-shifted light is propor¬ 
tional to the variance of the checkerboard density fluctu¬ 
ation ( 0 ^) — ( 0 )^, and thus to the static structure factor 
5'(kcb) = / dw ^(kcb, w). We show the variance of den¬ 
sity fluctuations as fllled symbols in Fig. 3, and observe a 
divergence when approaching the critical point P/Pa = 1 
from either side, heralding a second-order phase transi¬ 
tion. The inset displays the variance of the density fluctu¬ 
ations on a double logarithmic plot to illustrate the scal¬ 
ing behavior. The variance is plotted as a function of the 
Hamiltonian coupling parameter A, derived from mea¬ 
sured quantities and using a theoretical model (see Sup¬ 
plementary Note 1). In the normal and the self-organized 
phase we extract critical exponents of 0.7(1) and 1.1(1), 
respectively. We attribute the deviation from our previ¬ 
ous measurement in the normal phase, which gave 0.9(1) 
[ 2 ], to the refined model used for the scaling of the hor¬ 
izontal axis and the improved measurement scheme that 
allows us to directly distinguish between density fluctu¬ 
ations and a density modulation. Current theoretical 
research taking into account the open character of the 
system due to cavity dissipation predicts an exponent of 
1.0 [a [M]. The difference between the experimentally 
observed exponent and the predicted value might origi¬ 
nate from finite size effects and the presence of a small 
symmetry breaking field. Further, the theory models do 
not include damping of the momentum excitation due to 
atom-atom collisions [551l37| . 

Characterization of the quasi-particle mode 

The access to the dynamic structure factor S'(kcb,a;) 
allows us to characterize the quasi-particle mode which 
emerges due to the long-range interactions in the gas. 
When adiabatically switching on the long-range interac¬ 
tions (P ^ 0 ), new quasi-particle modes of polaritonic 
character form, where intra-cavity photons are admixed 
to the recoil momentum states. A diagonalization of the 
Hamiltonian for P ^ 0 leads to the definition of a quasi¬ 
particle mode for the interacting system with annihila¬ 
tion and creation operators c and c^, respectively (Sup¬ 
plementary Note 1). From our measurements. Fig. 2, we 
can directly extract the energy Hug of this quasi-particle 
mode as a function of P/Pcr- To this end, we fit a reso¬ 
nance curve of a damped harmonic oscillator to both side¬ 
bands of ^(kcbja;) (Supplementary Note 2), whose peak 
positions corresponds to Wg, see Fig. 4. We observe the 
mode softening towards the critical point from both sides 
of the phase transition. The width 7 of the sidebands is 
displayed in the lower panel of Fig. 4 and characterizes 
the damping of the quasi-particle mode. For our param¬ 
eters, the main constituent of the quasi-particle mode 
is the atomic component, while the light field is only 
weakly admixed. We thus attribute the observed damp- 
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Figure 4. Characterization of the quasi-particle mode. 

Frequency cJa (a) and decay rate 7 (b) of the quasi-particles as 
a function of relative pump power P/Pcr, as extracted from 
the fit to the dynamic structure factor, Fig. 2. The grey 
shaded area in the top panel results from an ab initio calcula¬ 
tion of the expected soft mode frequency, taking into account 
the experimental uncertainties in the determination of the 
coherent cavity field and the depth of the optical lattice re¬ 
sulting from the transverse pump field. Close to P/Pcr = 1, 
for Ws/(S tt) < 400 Hz, the uncertainty in modeling atom loss 
leads to a substructure which we omit in the graph. The solid 
line in the lower panel is a fit with a phenomenological func¬ 
tion to the data (see Supplementary Note 3). Vertical and 
horizontal errorbars indicate the statistical errors (s. d.) re¬ 
ported from the fit, and the error (s.d.) in the determination 
of the critical point, respectively. 

ing mainly to the decay of atomic momentum excitations. 
The finite decay rate n of the cavity light field gives rise 
to an additional damping of the quasi-particle mode es¬ 
timated to be only a few Hertz (Supplementary Note 1). 
The behavior of the damping rate 7 has been studied the¬ 
oretically and originates from a resonant enhancement of 
the Beliaev damping of the checkerboard density-wave 
[55fl37| . Our characterization of the quasi-particle mode 
is consistent with earlier measurements mm, but now 
also extends into the organized phase because we can dis¬ 
tinguish density fluctuations from a density modulation. 


Occupation of the quasi-particle mode 

The occupation (c^c) of the quasi-particle mode can 
be extracted from the observed sideband asymmetry. As 
can be seen in Fig. 2, the red-shifted sideband dom¬ 
inates over the blue-shifted one. For a system in its 












5 



Figure 5. Number of quasi-particles. Number of quasi¬ 
particles (c^c) — nr as a function of relative pump power 
P/Pcr, extracted from the sideband asymmetry in the dy¬ 
namic structure factor. The grey shaded area shows the result 
from an ab initio calculation of the expected quasi-particle 
mode occupation (Supplementary Notes 1 and 3), taking into 
account the experimental uncertainties in the determination 
of the coherent cavity field and the depth of the optical lat¬ 
tice resulting from the transverse pump field. Shown as black 
dashed line is the calculated thermal occupation nr of the 
quasi-particle mode due to the finite temperature of the BEC 
for a temperature of 38 nK. Vertical and horizontal errorbars 
indicate the statistical error (s. d.) reported from the fit, and 
the error in the determination of the critical point (s. d.), 
respectively. The strongly increased vertical errorbars close 
to P/Pcr = 1 arise from the decreasing sideband asymmetry, 
while their individual errors stay roughly constant. 


ground state, only creation processes are possible, lead¬ 
ing to a vanishing blue-shifted sideband. This has been 
used for thermometry of trapped ions and cavity optome¬ 
chanical systems |3H1 [31] ■ From the observation of the 
finite blue-shifted sideband in our experiment we infer 
that the system is in a steady state close to its ground 
state. The continuous measurement process via cavity 
decay constantly creates and annihilates quasi-particles 
at rates 2K{6a^Sa)- and 2K{Sa'^6a) + . Here, {Sa^Sa)± is 
the integrated spectral weight of the blue (-I-) and red (—) 
sideband, respectively. This measurement backaction ef¬ 
fectively gives rise to a heating rate of the system, and 
thus to a finite occupation of the quasi-particle mode. 
At the same time, the finite decay rate 7 of the quasi¬ 
particle mode changes this occupation: On one hand, 
quasi-particles will be annihilated due to this dissipa¬ 
tion channel at rate 2j{c^c}. On the other hand, it cou¬ 
ples the quasi-particle mode to a thermal heat bath pro¬ 
vided by the atomic cloud. This creates quasiparticles 
at rate 2^nT, where nx is the thermal occupation of the 
quasi-particle mode calculated from the Bose distribution 
function, and T = 38(10) nK is the temperature of the 
BEC, measured independently from absorption images. 
In steady state, the different contributions are balanced 
according to the rate equation (see Supplementary Note 


1 and Supplementary Figure 1), 

2k ((da^(5a)_ — {5a)5a)+) = 27 ((c^c) — nx) ■ ( 2 ) 

We can determine the occupation (c^c) — fix of the 
quasi-particle mode, using the weights of the sidebands 
and the dissipation rates 7 (empirical fit) and k (see 
Fig. 5). We observe an average occupation of the quasi¬ 
particle mode on the level of only a few quanta. In the 
organized phase an increase of the mode occupation to¬ 
wards the critical point seems visible. The occupation is 
expected to diverge when approaching the critical point 
since the energy of the soft mode vanishes and the atomic 
damping rate goes to zero IMIEZI- This situation is very 
similar to the enhanced thermal occupation of roton-like 
states predicted for dilute quantum gases with dipolar 
interactions at finite temperatures [9|. 

DISCUSSION 

We used vacuum-stimulated scattering of light to di¬ 
rectly measure the dynamic structure factor of a quan¬ 
tum gas with cavity-mediated long-range interactions. 
Access to the dynamic structure factor allowed us to 
characterize the relevant quasi-particle mode while the 
system crossed a structural phase transition, to distin¬ 
guish density modulation and density fluctuations, and 
to measure the critical exponents of the density fluctua¬ 
tions. We further extracted the finite occupancy of the 
quasi-particle mode under the influence of measurement 
backaction due to cavity decay and an atomic bath at 
finite temperature. While this measurement was moti¬ 
vated by the very specific setup used to create cavity- 
mediated long-range interactions, an extension to more 
general settings seems possible. The approach of apply¬ 
ing quantum optical methods based on strong matter- 
light interaction to the investigation of dilute ultracold 
gases offers unique possibilities for the non-destructive 
real-time investigation of quantum matter and its phase 
transitions 

METHODS 

After centering an almost pure ®^Rb Bose-Einstein 
condensate (BEC) trapped in a crossed-beam dipole trap 
with respect to the TEMqo cavity mode, the transverse 
pump power P, at wavelength Ap = 785.3 nm, is in¬ 
creased over 100 ms to a relative coupling strength of 
P/Pcr « 0.46. Subsequently, the power P is linearly 
increased over 0.5 s to P/Pcr ~ 1.38, while the stream 
of photons leaking out of the cavity is detected in a 
balanced heterodyne configuration using a local oscil¬ 
lator (LO) power of 2.2 mW and balanced photodiodes 
(Thorlabs PDBllOA) [3D|. The extracted quadratures 
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at a beat frequency of 59.55 MHz are mixed down to 
50 kHz, amplified, low-pass filtered and digitized using 
high-speed analog-to-digital converters with 2 us resolu¬ 
tion (National Instruments PCI-6132). The response of 
the heterodyne system is 2.2 per cavity photon. 

The temperature of the initially prepared BEC was de¬ 
termined from absorption images to be T = 20(10) nK. 
The far-off resonant transverse pump beam heats the 
BEC during probing to a temperature of T = 38(10) nK 
at the critical point. Residual atom loss of 26% during 
probing is included by rescaling the relative coupling axis 
according to the proportionality Per oc N~^. 

The phase transition point is characterized by a steep 
increase of the intracavity photon field. We fit the rise of 
the photon field once it has first exceeded a mean intra¬ 
cavity photon number of 4.5 with a saturation function 
Po ■ {1 — tcr/ty^ ■ With the extracted occurence time of 
the phase transition, we can convert the time axis into a 
relative coupling axis P/Pa. The relative statistical er¬ 
ror of Per according to this procedure is given by 5 • 10“^, 
including intensity fluctuations of the transverse pump. 
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SUPPLEMENTARY INFORMATION 


Supplementary Note 1: Theoretical description 


Relation between the structure factor and the power spectral density in our system 


We start with the definition of the dynamic structure factor as the spatial and temporal Fourier transform of the 
correlations of the density fluctuations 6p{r, t) [ST] 




1 


2t:N/V , 


drdt {5p{v, t)5p{Q, 0)), 


where V is the volume. With the relation p(r) = y this can be rewritten in Fourier space as 


5(k,cc) = 


1 


27riV 


dte*-‘(<5pk(t)(5p_k(0)). 


(3) 


(4) 


In steady state, the cavity light field operator a is given (neglecting small variations in the dispersive shift) by [S2 


a = 

Ac — in 

where the order parameter 0 = 0o + <50 

0 = 00 + [ d^rcoskxcoskzSp{Y,t) 


— 00 + ^ {^P{k,0,k) + Sp(^k,0-k) + Sp(^_k^o,k) + Sp(^_kfi-k)) ■ 


(5) 

( 6 ) 
(7) 


is written as a sum of a static density modulation 0o = (0) and time-dependent fluctuations. Here, the indices 
indicate the momentum along the x-, y-, respectively z-axis. We choose the coordinate system such that the cavity 
mode is oriented along the x-direction, and the transverse pump beam propagates along the z-direction. 

We now expand the cavity light field around its mean field solution a, a = a + So, (and correspondingly the order 
parameter as 0 = 0o -f <50), and use equation to find for the correlator 


{d\t)a(0)) = |ap + (M^(t)M(0)) 


A? 


00 + YQ{dp(k,0,k)5p{-k,0-k) + Sp(^kp-k)dp(-k,0,k) 


+dp(-k,0,k)dp(k,0,-k) + Sp(^-k,0-k)dp{k,0,k)) 


( 8 ) 


A? 


®0 + ^{^P{k,0,k)dp(-k,0-k)) 


Here, we made use of the fact that {6pq{t)dpq/{0)) = 0, unless q = —q', and that the system is symmetric under the 
transformations x —> —x and z —>■ —z due to the involved standing waves. Since the Fourier transform of the temporal 
correlation function (a^(t)a(O)) is the power spectral density of the intra-cavity light field, we finally obtain a direct 
relation between the dynamic structure factor and the spectrum of the intra-cavity light field: 


PSD{lj) = J dte-*‘"‘(a^(t)a(0)) 

n 2 




27r0Q(5(a;) + ^^S'(kcb,a;) 


(9) 

( 10 ) 


The power spectral density PSD{uj) thus quantifies both the checkerboard density modulation (first term, zero 
frequency bin) and the dynamic structure factor at the wavevector kcb = (±fc, 0, ±fc) (second term). 
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Finally we can extract the dynamic structure factor as 

5(k,b,a;) = ■ 


( 11 ) 


For the experimentally relevant case of two nearly degenerate, circularly polarized cavity modes, we find accordingly 

1 


5(kcb,w) = 




rjf + r ]2 N \2Tr 


-PSD{uj) - |ar(5(w) 


( 12 ) 


where the definitions of rji and Ac^eff are given below, and a and PSD{uj) refer to the according superposition of the 
bare cavity modes coupling to the atoms. 


Derivation of the Hamiltonian describing the fluctuations of the system 


In this section we derive an effective Hamiltonian describing the fluctuations of the system starting from the many- 
body Hamiltonian, 


l^mb — He + i?a + l^a-c + HgB > 


(13) 


with 


He 

H. 


2 

-hAc'^ala^ 


d^r'Fl(r) 


r r ^2 


2 m 


+ Vp cos^{kz) + ^'kl(r)'F(r) 


^(r) 


i?a 


= ^ / d^rd'l(r) hrji cos{kx) cos{kz){ai + a-) + HUq cos^{kx)a\ai 'l'(r) 


2 

Hsb = K X! ^*(®* + ) ■ 
2=1 


(14) 


Here we take into account that the cavity is supporting two degenerate (neglecting a small birefringence on the order 
of the cavity line width), circularly polarized cavity modes with annihilation operators ( 01 , 02 ), which are coupled to 
the atoms with two-photon Rabi frequencies ( 771 , 772 ) |S3j . The maximum dispersive shift of the two cavity modes due 
to the dispersive coupling of a single atom is described by (C/q , C/g). A symmetry breaking term ilsB is introduced, 
which is proportional to a real-valued effective cavity drive amplitude C, (see also definitions in |S2jl. 

Following references [S31IS4| . we expand the atomic and cavity field operators (4/, oi) around their mean-field values 
(V'o^ao), 


T = (VAt/^o -F (15) 

ai = ap + dcLi. (16) 


Expanding Hamiltonian Flmb, Eq. (131, in the limit |Ac| ^ k up to second order in the fluctuation operators 
(ildt, Mi), we find the quadratic Hamiltonian 


H-(2) = iJo + 71^ (M, + MJ) + , 


(17) 
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with 


Hr, = / d^r 


2m 


{dl + dl) + V^{z) + z)(a^ + a^*) 


2=1 


+^Xl^o(a;)|ao 


<5^' + ^ 52D V'o + 2g2B\il’o\^S¥6i>, 


(18) 


where A* = Ac — BqUq. We introduced here the spatially dependent classical lattice potential Vp{z) = VpCOS^[kz) 

of the transverse pump, the spatially dependent two-photon Rabi frequencies r]i{x,z) = cos{kx) cos(fcz) with 

maximum pump Rabi frequency ftp, detuning A^ = ujp — uja between atomic resonance frequency and pump light 

frequency, and single-atom coupling strengths Further, the light shift per photon is Ul{x) = cos^{kx) . Along 
the third direction, y, we assume a homogeneous system and use the according contact interaction strength g 2 D |S5| . 
We introduced the definition S& — y/N f d^r((5^'^ -I- cos{kx) cos{kz)ipo of the fluctuations of the order parameter 
0 = f cos(kx) cos(kz) around the mean field value ©o = / d^r\ipQ\'^ cos{kx) cos{kz). In a similar way, we 

define the fluctuation operator SB = yfN J d^r((5'I'l -|-(5'l') cos(kx)‘^ipo of the bunching operator B = f d^rd'l'l' cos(fca;)^ 
around its mean-field value Bq = J d^rlV'oP cos(fcx)^. 


The mean-field values are determined by the corresponding Gross-Pitaevskii equation derived from H^^^: 






+ dl) + Vp{z) -f U^o{x)\a. 


2m 

2 


i |2 

ol 


2=1 


+^^Yv^ix, z){al + al* -by2D|^/’on Mx,z), 


2=1 


(19) 


with Og = , assuming adiabatic following of the cavity field. 

To find the collective excitations of the system around the mean-field solution, we expand in Bogoliubov modes 
hj (which diagonalize Hq with eigenenergies Ej) with amplitudes Uj{r) and Vj{r), 

6it{r) = Y {uj{r)hj + v*{r)hfj , (20) 

3 


resulting in 


3 

2 ^ ^ 

+ kY^ -^ISalSdi. 

i=l 

Here, we introduced interaction matrix elements x] = {ipo\cos{kx) cos{kz) + cos"^{kx)\uj + Vj), describing the 

overlap between the ground state wave function and the Bogoliubov excitations as presented in reference [S3] . From 
numerical calculations we know that only a single Bogoliubov mode ko = k with energy Eq = E is dominantly 
contributing via the matrix elements Xg = x*) all other matrix elements are suppressed by more than two orders 
of magnitude [S3]. This leads to 


2 

Ej^kj + kY^Vi {hj yj ^Xj* + h-c.^ -|- Sdi^ 

i=l 


2 

=Ek'fk + kY -K 


2 

Ef^Y, 

i=l 



( 22 ) 


If we now use the definitions wq = A/S for the bare energy of the uncoupled system, \i = V^T7i\/|x*l for the 
coupling strength and 5b = keE* for the collective atomic fluctuations (using y/^ = \/|x*|e*'^’ « a/|x*|, where 
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(jii « 0, i.e. (j)i « (^ 2 ), we finally arrive at the fluctuation Hamiltonian 

2 2 

H = hoJoSPSb + Ai^5ai+(5aJ^ ^(5S + <56^j . (23) 

i=l i=l 

In order to further simplify this expression, we apply a transformation which reduces the Hamiltonian describing two 
cavity modes to a Hamiltonian involving coupling to only one effective cavity mode. We introduce two new cavity 
modes with fluctuation operators ( 61 , 02 ), which are a linear superposition of the original modes {dai,Sd 2 ) 

Sdi = Adi + Bd 2 (24) 

662 = C 61 + D 62 ■ (25) 

Applying this transformation and requesting that the new operators obey bosonic commutation relations ([ 61 , oj] = 

1 , [ 62 , 62 ] = 1 ), that mixed terms of the new operators vanish, and that mode 62 decouples from the atoms, we find 

the coefficients 


A = 


AiA 


( 2 ) 


A 2 A 


( 1 ) 


C, B = 




1 + 


A?A 


( 2 ) 


A?A 


( 1 ) 


c = - 




B = 


(it) 


^ +1 


1 + 


A?A)^) 

A?A<i> 


Inserting these coefficients into If yields 


(A? + A|)Ai^^^ q 


H = -n- 


AfA 


( 2 ) 


+ A|Ai^^ ^ ^ A^A 

hijXf + A^di + 6|)(66 + 6b^) + hujoSb^Sb. 


6|6i — h 


a?a(^)^ + a^a«^ 


,( 2 ) 


AiA 


( 1 ) 


-O 2 O 2 


(26) 


(27) 


The mode 62 is now decoupled from atomic motion and we can write the effective fluctuation Hamiltonian for the 


system 

H/h = ujo 6 b^ 6 b — Ac,eff6a^6a - 1 - Aefi((5a^ -|- 6 d)( 6 b^ + 6 b ), 

(28) 

where we used 

6 d = 61 

(29) 


A (Af- 1 -A 2 )Ac^^Ac^^ 

A?A(^)+AiAW 

(30) 


Aeff = a/Ai + A|. 

(31) 


Langevin description of the coupled system 


The goal of this section is to derive coupled quantum Langevin equations effectively describing our system. The 
hierarchy of the relevant energy scales determines their derivation. The fastest timescale in the system is given by the 
cavity field decay at rate k « 27r x 1.25 MHz, followed by the variable coupling rate Aeff, which reaches « 27r x 100 kHz 
at the critical point, and finally the damping rate 7 of the atomic excitation, which is on the order of a few hundred 
Hertz |S2| . Following this hierarchy, we first couple the damped optical cavity mode to the Bogoliubov mode 5b. This 
results in new polariton modes of the system which describe the quasi-particles of the long-range interacting system. 
Only thereafter we introduce the damping of the atomic polariton mode at rate 7 . This approach is justified as the 
admixture of the cavity mode to the polariton is very small (wq/Ac .eff ^ !)■ 


We start with Hamiltonian Eq. (28), 


H/h = ujoSb^6b + Ac,eff6a^6a -I- Aeff(<5a^ -I- 6d)(6b^ + 6b) 


(32) 
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The corresponding quantum Langevin equation with cavity held damping rate k read 


Sa = (—iAc,eff — K)Sa — iXesi^b + 6b) + (33) 

si) =—iuJoSi) — iXc.s{Sa + da') (34) 


The bosonic operator Uin in Eq. (33), describes vacuum input huctuations of the surrounding electromagnetic held 
modes which are characterized by the correlation functions {ain{t)al^{t')) = Sit — t') and {a\^(t)ain{t')) = 0. 

In order to hnd the new eigenmodes (polaritons) of the coupled system we rewrite the Langevin equations in matrix 
form, 


d 

dt 


/Sd\ 


Sd^ 

Sb 


= Mn 


/Sd\ 

Sd^ 

Sb 


+ 


/din\ 
t 


“in 

0 


with the matrix 


Mn = 


w) 

wJ 

V 

0/ 

/ ^^c,eff 

hi 0 

^Aeff 

^ Acff'^ 

0 

^^c,eff ^ 

^Aeff 

^Aeff 

Aeff 

zAeff 


0 

V ^Aeff 

^Acff 

0 



(35) 


(36) 


Diagonalization of Mq via the transformation 


S-^MnS = D = 


0 


V 


0 0 0 ^ 

^ b b 

0 —iuJs b 

0 0 iujsj 


with, neglecting terms of order 


S = 


1 

0 

zAeff 


0 

1 

^Aeff 

_ 


^ liAc,eff + «^ —2Ac,ef'f+«^ 


—-^eff 

^c,eff 
_ —Aeff 

Ac,eff + *« 
l-\-U}s/0J0 
2 

l — LJs/^Q 

2 


LOs 

^0 

(jJs 

UJQ 

UJs 

UJQ 


— Aeff 

Ac,eff — 

_ —Aeff 

Ac.eff+^Z* 

2 

2 


^0 

(jJs 

tJQ 

UJs 

UJQ 

UJs 

UJQ 


and 


5-1 = 


-Ae 


Ac,eff 

Aeff 

>Ac,eff-*« ' 


! ^0 
UJs 


Aeff 


Ac,eff+^« ^ 
_ — Aeff 

Ac,eff+^« ' 


/ UJQ 
UJs 

/ UJQ 


_ ^Aeff 

^Aceff+Z* 

iAeff 

—'iAc,eff+Kl 
l-\-UjQ/uJs I ijJs 

2 y ^0 

1 —ct;o /UJs I UJs 
2 V ^0 


^Aeff 


zAc,eff+'^ 

— ^Aeff 

— 2 Ac,eff+^ 

1 — UJq/UJs I LOs 
2 y txJQ 

l-\-UjQ/uJs j UJf. 
2 y UJQ 


allows us to dehne the polariton mode operators d and c. Here, 


(37) 


(38) 


(39) 


LUs — UJ^yJ\ (Aeff/Aeff, cr)^ ('^^) 

is the eigenfrequency of the quasi-particle mode we are interested in |S 6 | . Note that S~^ is normalized such that the 
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bosonic commutation relations for w and defined as 


fit 

c 


( 5d\ 
5a) 
5b 

Vfiv 




(41) 


are valid. This way we achieve the following expressions for the operators d (c), annihilating an excitation in the 
optical (atomic) polariton mode with eigenfrequency Ac,eff (ws): 


d = 5d-\- - — {5b + Sb^) 


/Wo . 

c = W -Aeff 

V Ws 

1 /wl 


—5a 




Ac^efF Ac^eff “t” J 

^1 + J^)Sb+{l-,l^)5p]. 


(42) 


Again, this is correct up to terms of order wo/Ac^eff- For Aeti —F 0 these equations yield d —)■ Aa and c —F A6, such 
that the polariton modes evolve into the bare modes for vanishing coupling. 


Langevin equations for the polariton modes 


We now can write down the Langevin equations for the polariton modes. Starting from equation (35) we use 


to find 


d 

dt 


f^\ 

Vcty 


= -^s- 

dt 


/5a\ 

Aat 

5b 

\5b^J 


= S-^MnSS- 


/5a\ 

' Aflt 

5b 

WJ 


/din\ 

VO/ 


(43) 


d = (—iAc^eff — K)d + 


KOii, 


c = —iuJsC + / —A, 


eff 


Ac,efF tK A^^eff “b IK 


(44) 

(45) 


The input noise terms in the equation for c correspond to the quantum backaction: although the cavity bath is at 
T = 0, the c-mode is driven by quantum noise from the open channel. 


Damping of the c-mode 

We now formally introduce a damping of the atomic polariton mode which effectively models its collisional in¬ 
teraction with the surrounding Bogoliubov modes of the atomic cloud. In principle, this mode is also damped by 
the cavity decay channel. However, since wo/Ac,eff ‘C 1, this damping rate is on the order of a few Hertz and will 
be neglected here. We define input noise operators with correlation functions {cm{t)c[-^{t')) = 5{t — t'){l -b ut) and 
(c|jj(f)cin(t')) = 5{t — t')nT- The thermal occupation number ut of the atomic polariton mode is modeled by the 
Bose distribution function, evaluated at the soft mode frequency tOs of the coupled system. Accordingly, the Langevin 
equation ([4^ becomes 
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C = {-ioJs - 7)c + — A, 


-Oh 




eff 


^c.eff ^c.eff H” , 


(46) 


Occupation number of the c-mode 


With this result, we can now calculate the expectation value (c’^(t)c(<)) for the number of quasi-particles. From the 
solution of Eq. (46), 


c(t) = [‘ (v^e,„ + 

u 0 7 V 7 ^c,efF ^ 


• ^^ 

IK Ac,eff + tK 


we hnd 


A 2 +k;2 


(47) 


(48) 


which diverges towards the critical point with an exponent 0.5 as a function of 1 — Aeff/Ac 


Spectrum of the light field 


The main observable in our experiment is the photon spectrum of the light held leaking out of the cavity. In order 
to hnd a relation between this spectrum and the c- and d-mode, we apply a backtransformation using the matrix S, 
and move to frequency space (for the dehnition of the Fourier transformation see [S2jl: 


da{uj) = d{uj) — Aeff 


Ac,eff 4“ tK /cUq 


Ac,eff + V 


7 ^ (c(w) -4c1'(-a;)) 


da'I' (w) = S (w) - Aeff (c’l' (w) -h c(-a;)) 


Ac,eff + K^\ 


(49) 

(50) 


Using the Langevin equation in Fourier space [S^, the expectation value of the cavity output photon spectrum can 
be written as 

(5at(^)5a(a;')) = {d^u})d{uj')} -7 - ((ct(aj)c(t^')) + (c(-w)c^(-^^')>) 

Ac,eft + K^^‘> ^ > 


\ 2 

^eff ^0 
Ac.efT + 


i 27<5(o; - J)nT + ^<5(c. - ..') 


V 


\iujs -t 7 — 


(51) 


27<5(o; - o.')(l + nr) + - ^') 


\iuJs - 1 - 7 - 1 - 


We set {d^ {uj)d{uj')) = 0, as the photonic polariton mode occupation is vanishingly small. In the last line of equation 


(51), the hrst term in the sum corresponds to the blue-shifted sideband, while the second term corresponds to the 
red-shifted sideband. The asymmetry in the amplitude of the sidebands is thus given by the different factors ftp and 
(1 -I- hr). The expectation value {5fd{ijj)5d{w')) corresponds to the experimentally observed power spectral density 
PSD{u}), excluding the coherent part at w = 0. 
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Supplementary Figure 1. Energy diagram of the driven BEC-cavity system. The bare states are denoted by |np) ® 
|uc,Ua), where rip is the photon number of the coherent transverse pump field (shown are the manifolds of states for two 
different rip), is the intra-cavity photon number and Ua is the number of atoms with (px,Pz) = ±hkx). The coupling 

which conserves parity of Ua + Uc softens the atomic excitation spectrum for increasing coupling strength Aefi. Admixtures of 
states to the coupled states are indicated. At the critical coupling strength Aeff, cr, the excitation gap huJs vanishes, triggering 
the quantum phase transition to a self-organized atomic state. The finite cavity decay makes this phase transition of non¬ 
equilibrium nature as the corresponding decay processes (wiggling lines) give rise to a depletion of the ground state of the 
Hamiltonian system. The exiting cavity photons at frequency ujp ± uJs carry real-time information about the dynamics of 
the quantum many-body system, corresponding to the creation (red line) and annihilation (blue line) of quasi-particles. The 
collective momentum excitations can decay at rate 7 , which results in an effective decrease of the annihilation processes via 
photon decay. This leads to the observed sideband asymmetry. 


Sideband asymmetry 


We can now explicitly calculate the sideband asymmetry in the cavity output photon spectrum, defined as the 
difference between the spectral weight of the red-shifted and the blue-shifted sideband. With the definition of the 
mean intracavity photon number in steady state [S7| 

-I pOC) poo 

{a^ a) = — / (a’^(a;)a(a;'))da;dw' (52) 

27 !" 7-00 7-00 

we call the integral over the blue-shifted sideband {6a^6d) + , and the integral over the red-shifted sideband {Sd^6d)-. 
The difference between the two sidebands, written as a rate of photons leaving the cavity is thus 
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2k (^{5d)5a)- — {5a^5a)+) = 

2tt J_^ A2 + k2 UJs (uJs + UJ)^ + 7^ 


^27(nT + 1) + 


2KXlg ojQ \ 


- r±.- .- ‘ 

StT J_^ A2 ^g + k2 Ws (Ws - w)2 + 72 


27717 ’ ■ 


2 KAgg 


Wo 


^Lff + 


^ r ^0 27 

277 J_^ ^ A 2 ^g + k2 Ws { oJs +07)2+72 


= 2k 


^0 ^eff 

A 2 ^g + «2 


(53) 


Comparing this result with equation (48 1 , we find a direct relation between the sideband asymmetry and the 
expectation value for the number of quasi-particles: 


2k ((Sd^Sd)- — (5d^5d}+) = 27 ((c^c) — ut) ■ (54) 

This equation can be interpreted as a rate equation: quasi-particles can be created either via thermal excitation 
at rate 27717 ’ or via the loss of a red-shifted photon out of the cavity at rate 2K{5dX5d)-. On the other hand, quasi¬ 
particles can be annihilated via the loss of a blue-shifted cavity-photon at rate 2K{5dX5d) + , or via damping of the 
c-mode at rate 27 (c^c). These different processes are illustrated in Supplementary Figured together with an energy 
diagram of the coupled system. 


Supplementary Note 2: Data evaluation 

The dynamic structure factor 

For the Fourier analysis, the detected signal is cut into subtraces of length 0.02P/Pcr corresponding to time 
windows of 11ms. For averaging purposes, each subtrace overlaps halfway with the neighboring ones. The two 
digitized heterodyne quadratures Qi and <52 are first electronically demodulated at 50 kHz. For each subtrace, we 
compute the Fourier transform of Qi -\-iQ2 using frequency bins of 90 Hz. We normalize the Fourier transforms to the 
number of data points in the subtrace and the frequency resolution, and average over all 147 experimental runs. The 
obtained spectral density 5(w) of the intracavity photons is converted into a power spectral density of the light field 
leaking out of the cavity via PSD = lOlog (5(w) • hcjXc ■ 2K/lmW) with the cavity decay rate k = 2Tr ■ 1.25MHz and 
the energy per cavity photon of hc/Xc- 


Sideband fitting procedure 


A symmetric low-frequency noise feature develops around the strong coherent field in the self-organized phase, which 
we attribute to technical phase noise between the signal and the LO beam. Sufficiently deep in the self-organized 
phase, its amplitude is proportional to the coherent field intensity. Before fitting, the scaled averaged technical noise 
is subtracted from the data in a region up to ±2 kHz. Then, the sidebands ap pea ring in the dynamic structure factor 


in Fig. 2 are fitted with two damped Lorentzian functions according to Eq. 


A- 
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^incoh 


(w) = BG 


((uj+UJ^)2+^2y 


((uj-uj )^+ 7 ^)^ ’ where BG is the heterodyne noise level with the signal path blocked, and A_|_/_ the sideband amplitudes. 

The imperfect modeling of the technical noise around the phase transition does not permit to start the fit of the 
sidebands at zero frequency. In order to be dominated by the sideband amplitude and not by noise, we thus restrict 
the sideband fitting routine to a finite frequency Vminfif This introduces a minimum excitation frequency that can 
be obtained through the fit. This effect can be observed in Fig. 4 for the quasi-particle frequency very close to the 
critical point: it is limited to 1 kHz. To account for the influence of the cut-off frequency iZmm.flt, we extract the fit 
results for from 810 Hz (900 Hz in the superradiant phase) to 1170 Hz in 90 Hz steps. Their mean is used as 
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the output of the fitting routine. To determine the error on the fit parameters we select the maximum of one of the 
three following error estimates: the mean individual fit error, the standard deviation of the fit results for the set of 
or the standard deviation of the fit results for keeping versus not keeping the quasi-particle frequency fixed 
to the value of the ab-initio theory. 


Supplementary Note 3: Comparison of ab-initio theory with experimental data 


The ab-initio calculations used in Fig. 1-5 are based on Hamiltonian Eq. (13) with no free fit parameter |S3| . They 


include the two degenerate cavity TEMqo modes with circular polarizations ei and £2 that are driven off-resonantly 
by the transverse laser beam, which is linearly polarized along the y-axis. For the atoms, prepared in the hyperfine 
state (E, tof) = with respect to a quantization axis pointing along the cavity axis (x-axis), the ratio of the 

corresponding two-photon Rabi frequencies is 771/772 = 2.66/0.82. Here, we denoted F as the total angular momentum 
and mp the magnetic quantum number. A single maximally coupled atom induces a maximum dispersive shift of the 
two cavity modes of Hq = 27r x 71 Hz and Hq = 27r x 22 Hz. 

We furthermore take into account the measured coherent intracavity photon number and the transverse pump 
potential that creates a lattice depth of 2.8Er at the critical point, calibrated using Raman-Nath diffraction |S8) . 
Systematic uncertainties of these quantities are estimated to be 5% each (see shaded regions in Fig. 4 and 5). The 
Gaussian envelopes of the pump and cavity fields along the transverse directions [S9] are effectively included by 
weighted averages of V/, 771^2 and over the spatial extent of the atomic cloud. Here, the Thomas-Fermi radii of 
the condensate in the external harmonic trapping potential is calculated in the presence of the changing transverse 
pump lattice and atom number, and is given by Ry, Rz) = (3.5, 8.3,6.9) pm at the critical pump strength, 
assuming an atom number of = 1.0 x 10®. 

We use our theoretical model to calculate the dynamic structure factor in Fig. 2 from the PSD using Eq. (1) in 
the main text with a conversion error of 5%, and to calculate the energy coq of the uncoupled system. The relative 
coupling strength A/Acr shown in the inset of Eig. 3 is calculated from the experimentally determined quasi-particle 
frequency Ws and the theoretically obtained frequency of the uncoupled system, coq, using A^/A^ = 1 — The 

energy of the quasi-particle mode shown in Fig. 4 is evaluated using Eq. (40). The theoretical number of quasi¬ 


particles in Fig. 5 is calculated using Eq. (53) and (54). The thermal occupancy of the quasi-particle mode shown is 

evaluated using np = ^exp[|^] — 1^ , where kp is the Boltzmann constant, and T = 38 nK. The atomic damping 

rate 7 is fitted with a phenomenological function based on a dispersion function: 7(0:) = Pq{x — \)/{p\ -I- (cc — 1)^)^^. 
The fit results in the normal phase are {po,Pi,P 2 ) = (—27r • 33.8Hz,0.039,1.02) and in the self-organized phase, 
{po,pi,P 2 ) = (27r • 299 Hz, 0.016, 0.649). This function is used in Fig. 5 to evaluate the theoretical expectation for the 
quasi-particle number (grey shaded area). 
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